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A packing of spherical caps on the surface of a sphere (that is, a spherical code) is 
called rigid or jammed if it is isolated within the space of packings. In other words, 
aside from applying a global isometry, the packing cannot be deformed. In this 
paper, we systematically study the rigidity of spherical codes, particularly kissing 
configurations. One surprise is that the kissing configuration of the Coxeter-Todd 
lattice is not jammed, despite being locally jammed (each individual cap is held in 
place if its neighbors are fixed); in this respect, the Coxeter-Todd lattice is analogous 
to the face-centered cubic lattice in three dimensions. By contrast, we find that many 
other packings have jammed kissing configurations, including the Barnes- Wall lattice 
and all of the best kissing configurations known in four through twelve dimensions. 
Jamming seems to become much less common for large kissing configurations in 
higher dimensions, and in particular it fails for the best kissing configurations known 
in 25 through 31 dimensions. Motivated by this phenomenon, we find new kissing 
configurations in these dimensions, which improve on the records set in 1982 by the 
laminated lattices. 

52C25; 52C17 



1 Introduction 

One of the key qualitative properties of a packing is whether it is jammed, that is, 
whether the particles are locked into place. Jamming is of obvious scientific importance 
if we are using the packing to model a granular material. Furthermore, it plays a central 
role in studying local optimality of packings, because one natural way to try to improve 
a packing is to deform it so as to open up more space. 

Jamming has been extensively studied for packings in Euclidean space. See, for 
example, Torquato and Stillinger [44] and the references cited therein. However, it has 
been less thoroughly investigated in other geometries. In this paper, we investigate 
jamming for sphere packings in spherical geometry, that is, packings of caps on the 
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surface of a sphere. Jamming has previously been studied for sphere packings in 
(see Tarnai and Caspar [41]), but there seems to have been little investigation in higher 
dimensions. 

A packing of congruent spherical caps on the unit sphere S"^ ^ in R" yields a spherical 
code (that is, a finite subset of 5"^' ) consisting of the centers of the caps. The minimal 
distance of such a code is the smallest angular separation between distinct points in the 
code. In other words, the cosine of the minimal distance is the greatest inner product 
between distinct points in the code. The packing radius is half the minimal distance, 
because spherical caps of this radius centered at the points of the code will not overlap, 
except tangentially. A spherical code is optimal if its minimal distance is as large as 
possible, given the dimension of the code and the number of points it contains. (Note 
that this notion of optimality is different from requiring that no more caps of the same 
size can be added without causing overlap. Neither of these two notions implies the 
other.) 

Spherical codes arise naturally in many parts of mathematics and science (see Cohn [9] 
for a more extensive discussion). For example, in they model pores in pollen grains 
or colloidal particles adsorbing to the surface of a droplet in a emulsion formed by two 
immiscible liquids. In higher dimensions, they can be used as error-correcting codes 
for a constant-power radio transmitter. Furthermore, many beautiful spherical codes 
arise in Lie theory, discrete geometry, or the study of the sporadic finite simple groups. 

A deformation of a spherical code is a continuous motion of the points such that the 
minimal distance never drops below its initial value. A deformation is an unjamming if 
it does not simply consist of applying global isometrics (that is, the pairwise distances 
do not all remain constant). A spherical code is called rigid or jammed if it has no 
unjamming. It is called locally jammed if no single point can be continuously moved 
while all the others are held fixed. 

For example, in the face-centered cubic packing of balls in , the kissing configuration 
(that is, points of tangency on a given ball) consists of the vertices of a cuboctahedron. 
This code is locally jammed, but it is not in fact jammed (see Conway and Sloane [19, 
page 29] or Proposition 3.5 below). However, it can be deformed into an optimal 
spherical code, namely the vertices of a regular icosahedron, and the icosahedron is 
then a rigid code with a higher minimal distance than that of the cuboctahedron. 

As this example shows, deforming a spherical code is one way to improve it. Some 
optimal codes are not jammed; for example, the best five-point codes in consist 
of two antipodal points and three points orthogonal to them, and the three points can 
move freely as long as they remain separated by at least an angle of 7r/2. Furthermore, 
computer experiments suggest that an optimal code can have rattlers, that is, points not 
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in contact with any other point, although no such case has ever been rigorously analyzed. 
However, despite these issues, rigidity is a powerful criterion for understanding when a 
code can be improved. 

Note that whether a configuration is jammed depends on the ambient space. For 
example, the vertices of a square are jammed in but not in . 

For infinite packings in Euclidean space, there are more subtle distinctions between 
different types of jamming (see Bezdek, Bezdek and Connelly [7] and Torquato and 
Stillinger [43]) based on what sorts of motions are allowed. For example, are all but 
finitely many particles held fixed? Are shearing motions allowed? However, these 
issues do not arise for packings in compact spaces. 

Nevertheless, jamming seems to be a more subtle phenomenon on spheres than it is in 
Euclidean space. In Euclidean space, there is an efficient algorithm to test for jamming 
(see Donev, Torquato, Stillinger and Connelly [21]) but on spheres we do not know 
such an algorithm. The difficulty is caused by curvature, which complicates certain 
arguments. For example, in Euclidean space every infinitesimal unjamming extends to 
an actual unjamming, as we will explain in Section 2, but the corresponding procedure 
does not work on spheres. 

Once we have developed the basic theory of rigidity for spherical codes, we will devote 
the rest of this paper to applying it to analyze specific codes. We will focus primarily 
on kissing configurations (that is, spherical codes with minimal angle at least 7r/3, or 
equivalently the points of tangency in Euclidean space packings), because they form a 
rich class of spherical codes and include many of the most noteworthy examples. An 
optimal kissing configuration is one with the largest possible size in its dimension. 

As mentioned above, the face-centered cubic kissing configuration is not rigid, but we 
will prove that all of the other best configurations known in up to twelve dimensions are 
rigid. Along the way, we will produce what may be the first exhaustive enumeration of 
these configurations in up to eight dimensions, as well as a complete list of the known 
examples in nine through twelve dimensions (although we suspect that more remain 
to be discovered). Above twelve dimensions, the calculations become increasingly 
difficult to do, even by computer, but we analyze certain cases that are susceptible 
to conceptual arguments. In particular, we show that the kissing configuration of the 
Coxeter-Todd lattice K\2 is not rigid, while that of the Barnes-Wall lattice Ai6 is, 
although both lattices are conjectured to be optimal sphere packings in their dimensions. 

We particularly focus our attention on 25 through 31 dimensions, because of the 
remarkably small increases in the record kissing numbers from each dimension to the 
next (see Table 3 in Section 7 for the old records). The best configurations previously 
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Dimension 


Kissing number 


Dimension 


Kissing number 


1 


2 


17 


5346 


2 




18 


7398 


3 


12 


19 


10668 


4 


24 


20 


17400 


5 


40 


21 


27720 


6 


72 


22 


49896 


7 


126 


23 


93150 


8 


240 


24 


196560 


9 


306 


25 


197040 


10 


500 


26 


198480 


11 


582 


27 


199912 


12 


840 


28 


204188 


13 


1154 


29 


207930 


14 


1606 


30 


219008 


15 


2564 


31 


230872 


16 


4320 


32 


276032 



Table 1 : The best lower bounds known for kissing numbers in up to thirty-two 
dimensions. Numbers in bold are known to be optimal (see Schiitte and van 
der Waerden [39], Levenstein [29], Odlyzko and Sloane [36] and Musin [34]). 

known were not even locally jammed, but we see no simple way to deform them so as 
to increase the kissing number. However, in Section 7 we show how to improve on the 
known records. We give a simple argument that shows how to beat them, as well as a 
more complicated construction that makes use of a computer search to optimize the 
resulting bounds. 

The new records are shown in Table 1. It is taken from Conway and Sloane [19, 
page xxi. Table 1.2(a)], with three exceptions: the entry for R^^ was out of date in 
that table (see [19, Chapter 5, Section 4.3]), the entries for M'^ and M^"* come from 
Zinoviev and Ericson [46], and the entries for M?^ through M?^ are new results in the 
present paper. See also Nebe and Sloane [35]. For the best upper bounds known in up 
to twenty-four dimensions, see IVIittelmann and Vallentin [32]. 

Figure 1 shows a logarithmic plot of the data from Table 1, normalized for comparison 
with 32 dimensions. One can see the local maxima corresponding to the remarkable 
fig, Barnes-Wall and Leech lattices in dimensions 8, 16 and 24, respectively. Note 
also that the growth rate of the known kissing numbers drops dramatically after 24 
dimensions. 



Qeometry & Topology, Volume 15 (2011) 



Rigidity of spherical codes 



2239 




Figure 1 : A plot of log tv, — 55 log T32 , where t„ denotes the current record 
kissing number in M" 



2 Infinitesimal jamming 



We know of no efficient way to test whether a given spherical code is jammed. In 
principle, it can be done by a finite calculation, at least if the points in the code have 
algebraic numbers as coordinates, by using quantifier elimination for the first-order 
theory of the real numbers (see Tarski [42]). (The proof relies on Roth and Whiteley [38, 
Proposition 3.2].) However, quantifier elimination is not practical in this case. 

On the other hand, there are much more efficient tests for a related concept called 
infinitesimal jamming, using linear programming (see Donev, Torquato, Stillinger and 
Connelly [21]). Given a code {xi, . . . ,xn} C S"^^ , imagine perturbing x,- to + eyt. 
Then 

\xi + eyi\^ = 1 + 2{xi,yi)e + 0(6^), 
where (•, •) denotes the inner product, and 

{xi + eyi,Xj + eyj) = {xi,Xj) + i{xi,yj) + {xj,yi))e + O(e^). 

Thus, to preserve all the constraints up to first order in e, we must have {xi,yi) = for 
all /, and {xi,yj) + {xj,yi) < whenever {xi,Xj) equals the maximal inner product 
in the code. An infinitesimal deformation of the code {xi , . . . , Xf^} is a collection 
of vectors y\, . . . ,yf^ satisfying these constraints. It is an infinitesimal rotation if 
there exists a skew-symmetric matrix <I> G W^" such that yi = <I>x; for all /, and a 
code is infinite simally jammed if every infinitesimal deformation is an infinitesimal 
rotation. (Recall that the skew-symmetric matrices are exactly those in the Lie algebra 
of SO{n).) Note that for an infinitesimal rotation, 

{xi,yj) = {xi,^Xj) = -{^Xi,Xj) = -{xj,yi). 
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Thus, an infinitesimal rotation does not ciiange any distances, up to first order. Tiie 
converse is false (consider a square on the equator in , with an infinitesimal de- 
formation moving two opposite corners up and the other two down), but it is true for 
full-dimensional codes: 

Lemma 2.1 Let y\, . . . ,yi^ he an infinitesimal deformation of a code {xi, . . . ,xi\i} 
in S"^^ such that xi,...,xn span M". If {xi,yj) + {xj,yi) = for all i and j, 
then the deformation is an infinitesimal rotation. 

Proof First, note that if a linear combination '^^i '^i^i vanishes, then a;j,- = as 
well, because 




= - ^ai{xi,yj) 
= -(^aiXi,yj 

^ i 

= -{0,yj)=0 

for all J (and only the zero vector is orthogonal to a set that spans M"). Thus, there is 
a well-defined linear map <I> such that ^Xi = . Furthermore, the identity 

{xi,<^Xj) = {xi,yj) = -{yi,Xj) = -{^Xi,Xj) 

implies that <5 is skew-symmetric, because it holds for a basis of M" and hence 
(m, $v) = v) for all M, V G M". □ 

Every infinitesimally jammed code is in fact jammed. This is not obvious: one cannot 
simply differentiate a purported unjamming motion to get an infinitesimal unjamming, 
without dealing with two technicalities, namely whether there is a differentiable unjam- 
ming and what happens if all the first-order derivatives vanish. However, it is true, as 
pointed out by Connelly [13, Remark 4.1] and by Roth and Whiteley [38, Theorem 5.7]: 

Theorem 2.2 (Connelly, Roth and Whiteley) Every infinitesimally jammed spherical 
code is jammed. 

The cited papers deal with the more general setting of tensegrity frameworks, in which 
movable points can be connected by bars (with fixed lengths), cables (with specified 
maximum lengths) or struts (with specified minimum lengths), and they prove that 
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infinitesimal jamming implies jamming in this setting. For the special case of spherical 
codes, we connect each point in the code to the origin using a bar, and we insert struts 
between neighboring points (that is, those at the minimal distance). 

We do not know whether every jammed spherical code that spans the ambient space 
is infinitesimally jammed. For tensegrity frameworks, the corresponding statement 
is not true: if we place bars along the edges of a regular octahedron, and use two 
additional bars to connect its center with a pair of opposite vertices, then the framework 
is rigid, but flexing the center orthogonally to the two adjacent bars is a nontrivial 
infinitesimal deformation. We have not found such an example for spherical codes, 
but we expect that there is one. By contrast, infinitesimal jamming is equivalent to 
jamming for periodic packings in Euclidean space (see Donev, Torquato, Stillinger and 
Connelly [21]). Specifically, if we perturb x to x + ey and x' to x' + ey' , then 

\{x + ey) - {x' + ey')\^ = \x - x'f + 2{x -x',y- y')e +\y- y'\^e^. 

The second-order term is always nonnegative, so nonnegativity of the first-order term 
suffices to produce an actual unjamming. (Deforming the underlying lattice complicates 
the analysis, but the result remains true; see [21, Appendix C].) What goes wrong in 
the spherical case is that + is no longer a unit vector and must be normaUzed, 
which causes the distances to decrease. 

This is not merely a technicality: there seems to be no simple method to turn an 
infinitesimal unjamming into an actual unjamming. Nevertheless, in all our examples, 
we have been able to accomplish this (with some effort). 

The linear programming algorithm for infinitesimal rigidity testing works as follows. 
By Lemma 2.1, to test whether a full-dimensional code is infinitesimally jammed, 
we need only check for each pair of points whether the distance between them can 
be changed. In other words, in an infinitesimal deformation, are the maximum and 
minimum of {xi,yj) + {xj,yi) zero for all / and j? For each / and j, this gives rise 
to two linear programming problems, because we are imposing linear constraints on 
the perturbation vectors yi, . . . ,yi\i and maximizing or minimizing the linear function 
(xi^yj) + {xj,yi). (Of course, when {xi,Xj) is maximal in the code, the definition 
of an infinitesimal deformation requires that {xi,yj) + {xj,yi) < 0, so maximizing 
this linear functional is trivial. However, the other cases are nontrivial.) The code is 
infinitesimally jammed if and only if the optima in all these linear programs are zero. If 
not, then solving the linear programs will produce an infinitesimal unjamming, provided 
that we also bound the coordinates of the perturbation vectors (to avoid unbounded 
linear programs). 

In some cases we are aided by symmetry, because we only need to check one representa- 
tive from each orbit of the action of the code's symmetry group on pairs of points in the 
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code. For example, if the symmetry group acts distance-transitively, then we only need 
to check one pair of points at each distance. Using the approach of Donev, Torquato, 
Stillinger and Connelly [21], we can even reduce to solving one linear program, at the 
cost of randomization. Specifically, consider maximizing the linear combination 



where the coefficients c, j are chosen randomly from the interval [—1,1]. With proba- 
bility 1 , this approach will produce an infinitesimal unjamming if one exists. Thus, if 
the optimum is zero, then we can be confident that the code is jammed, although this 
does not constitute a proof. 

For most of the examples in this paper, we give short conceptual proofs of jamming. 
However, for some cases we must rely on computer calculations. In these cases, we 
have given rigorous, computer-assisted proofs by using exact rational arithmetic via the 
QSopt_ex linear programming software of Applegate, Cook, Dash and Espinoza [1] 
and checking every pair of points in the code. 

3 The kissing configurations of root lattices 

We begin by proving that the root systems D„ (for n > 4) and E(,, E-j and are 
infinitesimally jammed, while A„ is not jammed (for « > 3). Note that these root 
systems are the kissing configurations of the corresponding root lattices. 

In this section, we will use spheres of diameter \/2 instead of 1, because that is 
standard for these root systems and makes the inner products integral. Note that the 
theory of infinitesimal jamming in no way depends on this normalization. 

The following elementary lemma will play a key role in the proofs: 

Lemma 3.1 Let C and V he spherical codes with CCD and with the same minimal 
distance. If C is infinitesimally jammed within the vector space it spans, then in any 
infinitesimal deformation ofV, the inner products between points in C are unchanged 
(up to first order). 

The interesting case is when C is lower dimensional than V . 

Proof Let x and y be points in C, and let u and v be their perturbations in an 
infinitesimal deformation of V . We write u = uc + u± and v = vc + v±, where uc 
and Vc are in the span of C while u± and v± are in the orthogonal complement of the 
span. 
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The orthogonal projections to the span of C yield an infinitesimal deformation of 
C . (Here we need C and V to have the same minimal distance, since otherwise the 
conditions on which inner products can increase will differ.) Thus, because C is 
infinitesimally jammed within its span, (x, vc) + {uc,y) = 0. Furthermore, (x, v^) = 
{y, u±) = 0. It follows that (x, v) + {u,y) = 0, as desired. □ 

Lemma 3.2 The A2 root system is infinitesimally jammed. 

The A2 root system is a regular hexagon, and it is easy to show that every regular 
polygon is infinitesimally jammed. This simple observation provides a useful tool for 
analyzing more elaborate configurations via Lemma 3.1. 

Proposition 3.3 The D4 root system is infinitesimally jammed. 

Proof The minimal vectors of D4 have norm 2 and the possible inner products 
between distinct minimal vectors are 0, ±1 and —2. First, note that the automorphism 
group of D4 acts transitively on pairs of minimal vectors with a given inner product, 
so without loss of generality we can consider just one pair of points at each distance. 
(This transitivity fails for D„ with n > 4, because there are two orbits for inner product 
0, but the tiiality symmetry of D4 collapses them to one orbit.) 

Furthermore, D4 contains A2, and by Lemma 3.1 the distances in a copy of A2 cannot 
change because A2 is infinitesimally jammed within its span. This takes care of all the 
cases except for a pair of orthogonal vectors. 

We now have to show that if (x, y) = then (x, y) does not change in any infinitesimal 
deformation. Again by the distance transitivity of the automorphism group, we may 
assume X = (1,1,0,0) and j = (1,-1,0,0). Let u = (1,0,1,0), v = (1,0,-1,0), 
w = (1,0,0, 1) and z = (1,0,0, —1) be other minimal vectors of D4. Denote the 
first order perturbations of x,y,u,v,w,z by x' ,y' ,u' ,v' ,w' ,z' ■ Note that x -\- y = 
u + V = w + z, and that (x, m) = (x, v) = {y,u) = {y,v) = 1. Thus, by the A2 
embedding argument, we have 

(x, u) + {x, u) = 0, {y, u) + (/ , u) = 0, 

(x,v') + (x',v) =0, (3;,v') + (y,v) =0. 

Adding these equations, we get 

(x + y,u' + V) + (x' + 3;', M + v) = 0, 
or (using x + y = u + v) 

{u + v,u + v) + (x' + y , X + y) = 0. 
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Since we know that {u, u') = 0, etc., we get (denoting the first order change {u, v') + 
(m',v) in (m,v) by 6{u,v)) 

6{u,v) + 6{x,y) = 0. 

Similarly, we have 

5{u,v) + 6(w,z) = 0, 
6iw,z) + 6ix,y) = 0. 

From these three equations, elementary algebra implies 6(x,y) = 5(u,v) = S(z,w) = 0. 
This completes the proof. □ 

The lengthy argument for the last case amounts to verifying that a square embedded 
within D4 cannot be infinitesimally deformed. Note that this cannot simply be settled 
using Lemma 3.1, although the square is indeed jammed within its span, because the 
minimal distance in the square differs from that in D4. If that argument worked, it 
would also prove infinitesimal jamming for D3 , which is not true. (The A3 and D3 root 
lattices are isomorphic to the face-centered cubic lattice, whose kissing configuration 
is not jammed.) 

Corollary 3.4 The D„ root system (for n > 4) and the E^, Ej and E% root systems 
are infinitesimally jammed. 

Proof These configurations have norm 2 and inner products 0, ±1 and ±2, the same 
as in D4 . We first deal with Ef, , Ej and E^ . Their automorphism groups act distance 
transitively, so it suffices to consider a single pair of points at each distance. The D4 
root system embeds in each of these configurations (in fact, its Dynkin diagram is a 
subdiagram), so without loss of generality we can assume the pair of points is in D4. 
Now combining Lemma 3.1 and Proposition 3.3 completes the proof. 

The same proof works for D„ with n > 4, with one exception, namely that there are two 
orbits of pairs of orthogonal vectors, so the group does not quite act distance transitively. 
Specifically, the stabilizer of (1, 1, 0, . . . , 0) cannot interchange (1, — 1, 0, . . . , 0) and 
(0, 0, 1, 1, 0, . . . , 0). However, in both cases, these vectors are contained in a copy of 
D4 (namely, the one in the first four coordinates), so we can complete the proof as 
before. □ 

The An root system is locally jammed, and for « = 2 it is in fact jammed, but the 
unjamming for « = 3 extends to higher dimensions. 

Proposition 3.5 For n > 3 , the root system is not jammed. 
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Proof We will demonstrate an explicit unjamming. We realize the A„ lattice as the 
set of integral vectors in the subspace {(xq, . . . ,x„) : = 0} of M""*"^, so each 

vector in the A„ root system has one 1 and one —1 among its coordinates. 

We begin with A3 , whose twelve vertices form a cuboctahedron with six square facets 
and eight triangular facets, and we choose a perfect matching using non-overlapping 
diagonals of the squares. In the left half of the diagram below, we have labeled matched 
vertices with the same label (from 1 to 6). Each matched pair v and v' satisfies 
(v, v') = 0, so we can perturb v to (v + £v')/\/T+? and v' to (v' + £v)/\/T+^, 
as shown on the right. 




To see that this is an actual unjamming, we must check what happens to vertices 
V and w satisfying (v, w) = 1. The matched vertices satisfy (v',w') = — 1 and 
(v',w) = — (v,w'), and hence 

(v + ev' w + ew') 1 — 

-5 '- = < 1 

as desired. 

For A„ with n > 3 , we simply perform this unjamming on the copy of A3 in the first 
four coordinates, while leaving all the other points unchanged. All we need to check is 

that the minimal distances between the moving points and the unchanged points do not 
decrease. Suppose v is in A3 , paired with v' as above, while w is outside A3 and 
satisfies (v, w) = 1 . Then v and w have a single nonzero coordinate in which they 
agree, and the other nonzero coordinate of w is not among the first four coordinates. 
Because (v, v') =0, the vectors v and v' have disjoint supports, from which it follows 
that (v',w) = 0. Thus, 

{v + ev',w) _ 1 ^ 
so we have a genuine unjamming of the A„ root system. □ 
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Note that the first-order changes in the inner products of 1 vanish, as they must 
according to Lemma 3.i (because each pair of vectors with inner product 1 spans a 
hexagon); instead, there is a second-order improvement. 

4 Kissing configurations in five, six and seven dimensions 

The root systems analyzed in the previous section are conjectured to be optimal kissing 
configurations up through eight dimensions. In fact, is known to be the unique 
optimal kissing configuration in (optimality was proved by Levenstem [29] and 
by Odlyzko and Sloane [36], while uniqueness was proved by Bannai and Sloane [2]). 
However, in five, six and seven dimensions, the , E(, and E-j root systems are not 
the unique optimal kissing configurations. This was first observed by Leech [27, 28]. 

In this section, we enumerate what we believe may be a complete list of optimal kissing 
configurations in these dimensions, up to isometry. We find two configurations in , 
four in , and four in MJ . Those in and , as well as two of the ones in , 
are all kissing configurations of previously known sphere packings. However, two of 
the configurations in M7 do not arise naturally from densest packings in MJ and appear 
to be new. 

We have verified by rigorous computation calculations that all these configurations are 
infinitesimally jammed. (It would very likely be possible to check this by hand, but it 
does not seem worth the effort.) To do so, we used version 2.6 of the QSopt.ex rational 
LP solver [1] to get exact solutions to the linear programs. Specifically, for each pair 
of points, we used the software to verify that their inner product can neither increase 
nor decrease in any infinitesimal deformation. Of course, it would be far more efficient 
to take into account the orbits of the automorphism group on pairs of points, but it was 
easier to let the computer check many cases than to compute and keep track of the orbits. 
The one tricky part of the calculation is that the coordinates of the points are not always 
rational numbers (and using fioating-point arithmetic would make the calculations 
unrigorous). Fortunately, there is an easy fix: by rescaling certain coordinates, in each 
case we can use rational coordinates and carry out all the calculations with respect to a 
nonstandard inner product defined by a rational matrix. In fact, we can take the matrix 
to be diagonal. In the supplementary information for this paper (see Appendix A), we 
provide explicit rational coordinates for all of the configurations studied in this section, 
together with the corresponding inner product matrices. 

4.1 Methods 

Conway and Sloane [18] describe a method that conjecturally generates all the best 
sphere packings in low dimensions. For dimensions from five to eight, it works as 
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follows. We start with the checkerboard lattice 

D4 = |x G Z'^ : ^x,- = (mod2)|. 

!=1 J 

It leads to a sphere packing in with balls of radius \/2/2, which is conjecturally 
the densest sphere packing in four dimensions, as well as the unique densest periodic 
packing. The discriminant group DI/D4 = (Z/2Z)^ has four elements, the zero 
class a = (0, 0, 0, 0) and three nonzero classes, which we represent by the deep holes 
b = (1,0,0,0), c = (1/2,1/2,1/2,1/2) and d = (-1/2,1/2,1/2,1/2). Conway 
and Sloane conjecture that the densest packings in dimension n (with 5 < n < 8) are 
obtained by fibering over D4, that is, by positioning translated copies of D4 over 
the densest packings in dimension n — 4, scaled so that their minimal distance is 
1. Each translation vector will be one of a, b, c or d, and adjacent points in the 
(«— 4) -dimensional packing will be assigned different translation vectors. This idea of 
fibering is a generalization of the construction of laminated lattices (see Conway and 
Sloane [15]). See also Cohn and Kumar [1 1] for other calculations using the methods 
of [18]. 

4.2 Dimension 5 : 40 points on S"^ 

In dimension 5, the densest packings are conjecturally obtained by stacking layers 
of D4 on top of each other, that is, by arranging translated copies of D4 along the 
integers Z. At each point of Z, one has to make a choice of which translate of D4 to 
use. In other words, such a packing corresponds to a coloring of Z by four colors or 
labels a, b, c, d, so that no two adjacent integers have the same color 

For instance, the D5 lattice is obtained from the following coloring: 




Because the affine symmetry group of acts as the full symmetric group on the four 
classes in D^/D/i, it does not matter which two distinct classes are used as labels here. 

The best kissing configurations known in dimension 5 are obtained from the local 
versions of these packings. The D5 root system comes from the following diagram: 
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The subscripts mean that we take the 24 minimal lattice vectors in the central D4 layer, 
and 8 translated lattice vectors in each adjacent layer (specifically, those surrounding a 
deep hole). We call this configuration Cj^oa- 

Of course, there is no need to choose the same translation vector on both sides of the 
central layer, and we can form a competing configuration Cs^oi as follows: 




This is the kissing configuration of the other three uniform 5 -dimensional packings 
described in [18]. It is genuinely different from C5^o«> because it lacks antipodal 
symmetry. Its symmetry group has size 384, compared to 3840 for C^^Qa- 

4.3 Dimension 6: 72 points on 

Next, we consider kissing configurations in six dimensions. The best packings known in 
are obtained by fibering over D4 using a hexagonal arrangement of translates in the 
plane. There are four uniform packings, described concisely by the colors associated 
to the hexagon of six translates around a central copy of D4 . For more details we refer 
the reader to Conway and Sloane [18], and Cohn and Kumar [11]. 

These packings have four distinct kissing configurations, shown in Figure 2, and a 
simple case analysis shows that every possibility in this framework is equivalent to one 
of them (up to symmetries of the hexagon and permutations of the labels). They are 
genuinely different, as their automorphism groups have sizes 103680, 3840, 2304 and 
384, respectively. 

4.4 Dimension 7: 126 points on 

The best packings known in are obtained by fibering over D4 in an optimal 3- 
dimensional arrangement, namely one of the Barlow packings. The Barlow packings 
are themselves obtained by stacking translates of the hexagonal lattice. Here, we will 
be concerned with the kissing configurations, which involve only three hexagonal 
layers, so although there are uncountably many Barlow packings (corresponding to 
three-colorings of Z), we will only need to consider two configurations: the hexagonal 
close-packing, in which the top and bottom layers are mirror images of each other, and 
the face-centered cubic packing, in which they are point reflections of each other. 

First, we consider the face-centered cubic arrangement. Without loss of generality 
we label the central sphere at the origin a. Once we have specified the labels in the 
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C6-72C C^-lld 

Figure 2: Kissing configurations in M.^ based on the hexagonal lattice 

central hexagonal layer, the labels on the minimal vectors in the adjacent layers will be 
completely determined (since they will each already have three neighbors with different 
labels). It is not hard to check that in the central layer, the hexagon surrounding the 
central copy of D4 must be colored either bcdbcd or bcbdcd, up to a permutation of 
the colors, because no other choices will extend consistently to the neighboring layers. 
This determines all the labels in the two diagrams shown in Figure 3, except for the six 
points labeled a[ . For those points, one can check that if their labels agree with the 
label assigned to the central sphere, then they add one to the kissing number; otherwise 
they add zero. Thus, they should all be labeled a to get a kissing number of 126. 

Accordingly, we get two local configurations of tight packings, as shown in Figure 3. 
The code Cj-uea is the £7 root system, but Ci-u6b is not the kissing configuration 
of any of the tight packings in MJ described by Conway and Sloane [18]. 

The two other conjecturally optimal kissing arrangements in come from the hexag- 
onal close-packing via a similar argument and are shown in Figure 4. The code C7_i26c 
is the kissing configuration of the packing Aj described in [18], while Cj^ued does 
not occur in any of the tight packings analyzed in that paper. 

These four spherical codes can again be distinguished by the sizes of their automorphism 
groups, which are 2903040, 46080, 103680 and 3840, respectively. 
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Figure 3: Kissing configurations in MJ based on the face-centered cubic packing 

5 Nine through twelve dimensions 

Up through eight dimensions, we are quite confident that the known kissing configura- 
tions are optimal and reasonably confident that the list we have provided is complete. 
However, the situation in nine or more dimensions is very different. For example, 
in the highest kissing number achieved by any lattice is 272 (by a theorem of 
Watson [45]), while the best kissing number known is 306. Furthermore, there is little 
reason to believe that the usual constructions based on error-correcting codes (see 
Conway and Sloane [19, Chapter 5]) will prove adequate in high dimensions. Even in 
M^", computer searches have led to intriguing new kissing configurations (see Elser 
and Gravel [24]), although so far they have not improved on the known records. 

The laminated lattice Ag achieves kissing number 272 , and it is the only lattice that 
does. Here, we prove two properties of this kissing configuration: it is not locally 
jammed, but it is the largest possible kissing configuration in that contains the E% 
root system as a cross section. 

To describe this code, we first recall the structure of the root lattice , which is used 
to produce Ag . It the union of two translates of Dg , namely 

£8 =D8U (D8+(l/2,...,l/2)). 
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C7-I26C C7_i26rf 

Figure 4: Kissing configurations in based on the hexagonal close-packing 

The minimal vectors of are the 112 vectors of the form zte, it ej with / 7^ j 
(where e,- is the ith standard unit vector), and the kissing configuration consists of 
these vectors as well as the 128 vectors of the form (±1/2, . . . , ±1 /2) with an even 
number of minus signs. 

The vector v = (1, 0, 0, . . . , 0) is a deep hole of E^, and the lattice A9 is generated 
by Es, X {0} and v x {1} in M^. Its kissing configuration consists of 272 vectors, and 
is the union of the kissing configuration S of Dg (which contains 144 vectors) and 
the set T of 128 vectors of the form (±1/2,. ..,±1/2,0) with an even number of 
minus signs. 

By Corollary 3.4, the subconfiguration S is itself jammed, so it is futile to move those 
points of the code, but we will show that the points in T are not even locally jammed. 
The reason is that they have no neighbors outside of E^ x {0} , so they are free to move 
orthogonally to E^ . 

More formally, given (x, 0) G T with x G \ Dg , we move it to (x cos 6, \fl. sin Q) , 
where — 7r/2 < 9 < ir/l. For every point y ^ Eg with y ^ x,we have {x,y) < 1 and 
hence 

( (x cos 9, V2 sin 9),{y, 0)) = {x, y) cos 9 < \. 
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Therefore, the only points (xcos 0, \/2 sin ^) can come too close to are those in 
D9 \ , that is, those of the form ite,- it eg with 1 < / < 8 . We do indeed run into 
problems from these points when \6\ is too large, but not when \9\ is small. The 
relevant inner product is 

{{x cos , ±ei lb eg) = lb - cos 9 lb 

When \9\ is sufficiently small, the right side is clearly bounded by 1 . Specifically, we 
obtain a bound of 1 whenever 

l„l . def ^ _i \/5 - \/2 

\9\ < a = tan w 0.1247r. 

Thus, the point (x, 0) is not locally jammed. 

Furthermore, we can even simultaneously move all the points of T by moving (x, 0) to 
{xcos9x, V^&in9x), where —a < 9x < a and 9x has sign (—1)*^'/^, with being 
the number of negative coordinates of x. All we need to check is whether these points 
come too close to each other. The inner product between two of them is simply 

((jccos 9x, Vlsin 9x) , (jcos 9y, Vl sin 9y^)={x,y) (cos 9x){cos 0y)+2(sin6'.,:)(sin6'3,). 

If {x,y) = 1, then x and y have differing signs in two coordinates (that is, kx and 
ky differ by 2). In that case 9x and 9y have opposite signs, so 

(cos 9x)icos 9y) + 2{sin9x)(sin 9y) < {cos9x)icos9y) < 1. 

On the other hand, if {x,y) < 0, then the inner product is at most 2(sin^;c)(sin^y), 
which is again at most 1 (because \9x\ and l^^l are less than 7r/4). 

Proposition 5.1 No kissing configuration in ttiat contains the E% root system as a 
cross section can have more than 272 points. 

Proof Suppose C is a kissing configuration that contains the points (x, 0) for x £ Eg, 
with \x\'^ = 2. The remaining points must be of the form {y cos 9 , \/2 sin 9) with 
y eR^ satisfying \y\^ = 2 and -7r/2 <9 < 7r/2. 

Given any such point y, there exists a minimal vector x £ Eg such that {x,y) > \f2. 
This claim amounts to knowing the depth of the deep holes in the Eg root system, which 
correspond to the cross-polytope facets of its convex hull. (See, for example, Conway 
and Sloane [17].) It follows that 6* = is impossible, since (j,0) would come too 
close to the point (x^ 0). In fact, because the points (;k;, 0) and (ycos 9^ \/2 sin 9) in C 
must have inner product at most 1, we must have -v/2cos0 < 1. Thus, \ sin9\ > l/\/2. 
Without loss of generality, we will focus on the points with sin > 1/ V2. 
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Given two distinct such points (y cos 6, \fl. sin d) and (z cos sin ip) , their inner 
product is again at most 1 . Therefore, 

1 > (y, z) cos 9 cos 93 + 2 sin sin 99 > {y, z) cos 6* cos if + I, 

from which it follows that {y, z) < 0. (Note that cos 9 and cos 99 cannot vanish, since 
then the sine term would be too large.) 

There are no more than 16 vectors in 5^ for which all inner products between distinct 
vectors are nonpositive (see Boroczky [8, Theorem 6.2.1]). Thus, we have shown that 
C contains at most 16 additional vectors in each hemisphere, for a total of at most 
240 + 2 • 16 = 272, as desired. □ 

As a consequence, the best kissing configuration in cannot contain the best one 
in R*^ (namely, E^) as a cross section. Note also that in any 272-point kissing 
configuration containing Eg, equality holds throughout the proof of Proposition 5.1, 
so there are only finitely many such configurations, each consisting of ^8 with cross 
polytopes sitting above and below some of its deep holes. However, the relative position 
of these cross polytopes may vary. For example, in Ag one of the cross polytopes 
consists of the points ite,- — eg with 1 < / < 8 . If // is any 4x4 Hadamard matrix, 
we could replace this cross poly tope with the points itv, — eg, where vi, . . . ,V8 are 
the rows of the matrix 

'H/2 " 
H/2_ ' 

The resulting code is genuinely different, because it contains two points with inner 
product —3/2, while only ±2, ±1, ±1/2 and occur in the A9 configuration. 

Table 2 lists the best kissing numbers known in dimensions 9 through 12 (see Conway 
and Sloane [19, pages 139-140]). The configurations are constructed using constant 
weight binary codes, and in fact the E-j and Eg root systems can also be constructed 
in this way. Every binary code B of block length n, size N, constant weight 4 and 
minimal distance 4 yields a periodic packing in W , namely all the vectors in Z" that 
reduce to codewords in B modulo 2. The vectors of norm 4 in this packing form a 
kissing configuration of size 2n + 16N , consisting of the points ib2e; together with 
signed codewords from B (that is, vectors with arbitrary ±1 entries in the support of a 
codeword). 

For n = 9, 10, 11 and 12 one can achieve N = 18, 30, 35 and 51, respectively. 
These codes are unique up to isomorphism when « = 9 or 10, there are 11 of them 
for « = 11, and there are 17 of them for n = 12 (see [37, Table I]; the n = 9 case 
was proved by Ostergard [37], the n = 10 case by Barrau [3], and the n = II and 12 
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Dimension Best known kissing number Packing 

9 306 

10 500 Pm 

11 582 Pii, 

12 840 Pna 

Table 2: The best kissing numbers known in dimensions 9 through 12 

cases by Best [4, 5, 6]). Data files giving coordinates for these codes are available in 
the supplementary information (see Appendix A). 

For n > II there are multiple codes, and each of these binary codes yields a distinct 
kissing configuration. To see why, first observe that each contains a cross polytope 
consisting of the vectors of the form ±2e; . This cross polytope is uniquely distinguished 
by the property of having large valencies. Specifically, when ?i = 1 1 it consists of 
exactly the points with at least 88 neighbors (that is, points at inner product 2), and 
when n = 12 it consists of those with at least 136 neighbors. These distinguished 
cross polytopes must correspond under any isomorphism between two configurations, 
which must therefore be a signed permutation of the coordinates. However, it would 
then yield an isomorphism of the underlying binary codes. 

In nine through twelve dimensions, these kissing configurations are all jammed, due to 
the following proposition, whose hypotheses can speedily be checked by a computer 
calculation: 

Proposition 5.2 Let B be a constant weight code of block length n, weight 4 and 
minimal distance 4. If B has the following two properties, then the corresponding 
kissing configuration in W is infinitesimally jammed: 

(1) For i <j, there is a codeword x = (xi , . . . , x„) G B such that Xi = xj = I. 

(2) For every x ^ B and i such that x,- = 0, there exists y ^ B such that y, = 1 and 
d{x, y) = 4 (that is, the supports of x and y overlap in exactly two coordinates). 

Proof Let C be the kissing configuration obtained from B. First, observe that for 
/ 7^ j, the points 2e,- and lej are part of a scaled D4 root system embedded in C. 
Specifically, by (1), there exist k and i such that it ej ±€^±61 are all in C ; then 

{±2ei, ±2ej, ±2ek, ±2ee, ±ej ± ej ±ek± ei} 

is a scaled D4 root system. Therefore, by Lemma 3.1 and Proposition 3.3, the inner 
products in D4 and thus in the cross polytope {±2e, : 1 < / < cannot change in 
any infinitesimal deformation. Since the cross polytope spans M" , we can assume that 
the points in it are fixed. 
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Let 

X = Sid + SjCj + SkCk + s^eii £ C 

be any of the remaining vectors, with Si, sj, Sk, S£ E {±1} , and let y be the infinitesimal 
perturbation of x. Then x and let generate a scaled copy of the A2 root system in C, 
which is also infinitesimally jammed. Therefore the inner product between x and lei 
does not change to first order, so 3^,- = 0. Similarly yj = y^ = yi = 0. 

Now let p {i,j,k,£}. By (2), there exists b £ B with bp = \ and such that b 
overlaps the support of x in exactly two positions. Assume without loss of generality 
that b is supported in positions i, j, p and q. Then zte, it ej ± ep ± is in C for 
all combinations of signs. The four points x, z = —SjCi — sjej — Cp — Cq, Icp and 
Icq span a scaled copy of the D4 root system in C ; indeed, they form a D4 Dynkin 
diagram: 




2ep Icq 



Therefore, the inner product between x and Icp does not change to first order, so 
= 0. It follows that J = 0, and thus C is infinitesimally jammed. □ 

The situation is quite different for lattice kissing configurations. The best one known 
in is that of the lattice Aio; it contains 336 points and is not locally jammed (as 
one can check by a calculation like that for A9). The best lattice kissing arrangement 
known in M'^ comes from the lattice A™"; it contains 438 points and is also not 
locally jammed. 

The best lattice kissing configuration known in is 

that of the Coxeter-Todd lattice 
K\2, which is not a laminated lattice. Unlike what happens in the previous three 
dimensions, this 756-point code is locally jammed and has a transitive symmetry 
group. However, we will show that it is not in fact jammed. Aside from the A„ root 
systems, the Coxeter-Todd kissing configuration is the only example analyzed in this 
paper that is locally jammed but is not jammed. It is also remarkable because the 
Coxeter-Todd lattice is the densest sphere packing known in M^^ . This phenomenon of 
seemingly optimal packings with locally jammed yet unjammed kissing configurations 
also occurs in three dimensions, but we know of no other cases. 

To unjam the Coxeter-Todd kissing configuration, we will make use of its Eisenstein 
structure: K\2 is a Z[tj]-module of dimension 6, where is a primitive cube root of 
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unity. We view K12 as a subset of with the inner product {x,y) = 2Re(x,j)c, 
where {z, w)c = Ylj=i^j'^j the usual Hermitian inner product on C^. The minimal 
vectors of Ku have norm 4 with respect to (•,•). 

The key property of K12 is that for minimal vectors x and y satisfying x ^ ia/j, 
even the complex inner product (x,y)c is bounded in absolute value by 1. This is 
quite unusual and does not hold for most other Eisenstein lattices (or related types of 
lattices such as Hurwitz lattices), but it can be checked directly from the list of minimal 
vectors in Conway and Sloane [19, page 128]. More conceptually, it follows from the 
Eisenstein integrality of K\2 (that is, the fact that all the complex inner products are 
in Specifically, if |(x,y)c| > 1, then a quick enumeration shows that the only 

possibility is (x,j)c = ±ay(l — w), in which case (x, iba/j)c = 1 — w. However, a 
complex inner product of 1 — a; leads to a real inner product of 3 , which is impossible. 
Note that this property is an assertion about the minimal distance of the 126-point 
configuration in CP^ obtained by taking the quotient of this code by the action of the 
multiplicative group , which in this case just amounts to taking the quotient modulo 
the sixth roots of unity. In fact, it follows from Cohn and Kumar [10, Theorem 8.2] 
that the resulting code in CP^ is universally optimal. 

The Eisenstein structure breaks up the code into 126 hexagons (the orbits under 
multiplication by powers of the sixth root of unity — cj), and to unjam it we simply 
rotate each of these hexagons by arbitrary angles. If we rotate x to e^^x and y to 
e'^Py, then 

{^'x,e'^y)^ = e'''~^\x,y)^, 

and hence 

(/^x.e^'Py) = 2Re{e'^x,e''Py)^ < 2\{e'^x,e''Py)^\ = 2\{x,y)c\ < 2. 

If each hexagon is rotated by a slightly different angle, then this inequality will be strict 
whenever x and y lie in different hexagons. Therefore, we get an actual unjamming. 
Within each hexagon, the minimal distance has not changed. However, we can follow 
up such a deformation by moving the elements in each of the hexagons away from 
each other. For instance, decompose the hexagon into two equilateral triangles and 
move them away from each other, in a direction orthogonal to the plane of the hexagon. 
This process increases the minimal distance within each hexagon and results in a 
configuration of 756 points with minimal angle strictly larger than tt/3. Thus, we 
can not just unjam the Ku kissing configuration, but unjam it so thoroughly that no 
contact remains (much like in three dimensions). 

Note that it is difficult to predict results such as this based on dimension counting. For 
example, for large codes in S"^^ , the isostatic condition suggests that about 2{n — 1) 
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contacts per particle are needed to ensure jamming (see Torquato and Stillinger [44, 
page 2641]), because that yields N(n — 1) constraints for an A/^ -particle code, which is 
the same as the number of degrees of freedom. (Strictly speaking, we should subtract 
dim 0{n) = (2) from the number of degrees of freedom, but that is negligible when 
N is large.) For sufficiently generic constraints, one might expect this bound to be 
sharp, but it is far from sharp for highly symmetrical configurations. For example, 
each particle in the Coxeter-Todd kissing configuration is in contact with 82 others, so 
this configuration has far more contacts than the isostatic condition requires, but it is 
nevertheless unjammed. The A9 kissing configuration also has plenty of contacts, but 
it is not even locally jammed. 

6 Kissing configuration of the Barnes- Wall lattice 

In this section, we show that the kissing configuration of the Barnes-Wall lattice is 
rigid. This lattice is the densest sphere packing known in M^^, and it has the highest 
known kissing number. 

We will consider the 4320 minimal vectors of A16 as points on the sphere of radius 2, 
which is their usual normalization. Then the inner products between distinct elements 
of the kissing configuration lie in the set {±2, ±1,0, —4}. 

Recall that the minimal vectors of Aie are as follows (see Conway and Sloane [19, 
page 129]). There ai^e 480 of the form 2'/^(ibe,- ± ej) with / 7^ j, which span 
2^/^Di6, and 3840 of the form 2^'/^( ^^^^ ite,) , where / is the support of one of 
the 30 codewords of weight 8 in the first-order Reed-Muller code 7^A^(1, 4) of block 
length 16 and where the number of minus signs is even. Codewords in 7^A^(1,4) 
correspond to affine linear functions from F2 to F2 , with the codeword consisting of 
the values of the function at the 16 points in Fj. 

Proposition 6.1 The kissing configuration of the Barnes-Wall lattice A16 is infinites- 
imally jammed. 

Proof Consider any infinitesimal deformation of this code. The points 2^/'^(±ei ± 
ej) form a copy of the Dig root system (scaled by 2^/^ so the minimal distances 
match). This root system is infinitesimally jammed by Corollary 3.4, so after applying 
an infinitesimal rotation we can assume that all these points are fixed to first order. 
Furthermore, each pair of points with inner product ±2 spans a scaled copy of A2 , and 
each pair of antipodal points is contained in such a hexagon, so these inner products 
also do not change to first order. Thus, only the inner products and ±1 can possibly 
change. 



Qeometry & Topology, Volume 15 (2011) 



2258 



Henry Cohn, Yang Jiao, Abhinav Kumar and Salvatore Torquato 



We will now show that the facts in the previous paragraph together imply that the 
infinitesimal deformation must be identically zero. For consider a minimal vector 
X = 2~^/^ Y2i=i ^i^ki, where 1 < < • • • < < 16 and Si € {±1}, and suppose x 
is infinitesimally perturbed hy y = J2iyi^i- 

The inner product of x with the Dig minimal vector 2^^^{sieki + Sje^j) is 2 (note 
that the supports overlap and the signs cancel), and as noted above this inner product 
cannot change to first order. Because all the Dig vectors are fixed, we find that 
^iyki + sjykj = 0. This holds for every choice of ki and kj, and if ki, kj and ki are 
distinct, then the equations 

^iyki + Sjjkj = Sijk, + suki = Sjjkj + sukf = 

imply that yuj = yUj = Jki =0. Since this holds for all / = 1, . . . , 8, we see that the 
nonzero coordinates of x cannot change to first order. Next we have to show that the 
zero coordinates do not change either. 

Let the positions of the zeros be < • • • < ^8 > and consider the Di6 vector z = 
l^l'^iei^ + ^^„) for some m ^ n. It is orthogonal to x, and by Lemma 6.2 below, 
X and z are contained in a sublattice isometric to a scaled copy of D4 . Since D4 is 
infinitesimally jammed, the inner product (x,z) does not change to first order, and 
since Dig is fixed, z itself does not change to first order. Hence y^^ + yi„ = 0. 
Similarly, j^,,, - y^,^ = by replacing z with z! = l^^'^ifit^,^ - e^). Thus, j^?,,, = 
for all m, and x does not change to first order, as desired. □ 

All that remains is to show that every pair of orthogonal minimal vectors in A16 is con- 
tained in a scaled copy of D4 (Lemma 6.2). This would be trivial if the automorphism 
group acted transitively on these pairs, because it is easy to exhibit a pair for which 
it is true (from a copy of D4 inside Dig, for instance). However, there are actually 
two orbits. Each minimal vector in Aig is orthogonal to 1710 others, and under the 
stabilizer of that vector they form two orbits, of size 1680 and 30. This non-transitivity 
is actually quite remarkable, because it corresponds to a decomposition of the 4320 
minimal vectors into 135 cross polytopes: to form one of the cross polytopes, take any 
vector, its antipode and the stabilizer orbit of size 30. It is far from obvious that this 
works, but it forms a decomposition that is invariant under the automorphism group. 
By contrast, the 240 minimal vectors in form 15 cross polytopes, and the 196560 
in A24 form 4095 (see, for example, Elkies [23, page 6, footnote 3]), but in these 
more symmetrical cases there is no invariant decomposition. 

In principle, one could simply check Lemma 6.2 for the two orbits, but we will give a 
direct proof that requires less information about the action of the automorphism group. 
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Lemma 6.2 Any two orthogonal vectors of Aig are contained in a sublattice isometric 
toD4 scaled by 2^/^. 

Proof Let x and y be two orthogonal minimal vectors in Aig. Without loss of 
generality, we may assume that x = 2^/^{ei + €2), because the automorphism group 
acts transitively on the minimal vectors (see Griess [25, Theorem 10.2] and [26]). We 
will construct vectors w and z such that w,x,y,z form a D4 Dynkin diagram, up to 
scaling: 




In other words, {w,x) = {x,y) = {y,w) = and {w,z) = = (y,z) = —2. 

Then these vectors will span a copy of 2^/^D4. 

First, we analyze the case when y G I^/^Dk,. Then either y = ±2^/^{ei — €2) or 
y = 2^/^{±ei ± ej) with 2 < i <i.\iy = ±2^/^{e\ — €2), then we can assume the it 
sign is positive (by replacing y with —y if needed, since this does not change whether 
X and y are contained in a scaled copy of D4), and we can take z = 2^l^{e-i — e\) 
and w = 2'/^(e4 — £3). If y = 2'/^(ibe,- it ej) with 2 < i < j, then we can assume 
that y = 2^/^{ei ± ej) and take z = -2^^^{ei + ei) and w = 2'/2(ei - €2). 

In the remaining case, we have y = 2^^/^ Z]f=i ^i^ki with Si G {±1} and < k2 < 
■ ■ ■ < kg. Since x is orthogonal to y, either ki = I and ^2 = 2, or ^1 > 2. If k\ = 1 
and ^2 = 2, then ^'1+52 = 0, and without loss of generality we can take i'l = 1 and 
52 = — 1 ; then let w = 2^/^{s3ek3 — S4ek^) and z = —2^/^{ei + s^ek^). 

Finally, suppose ki > 2. This case is slightly more subtle, because z cannot be an 
element of 2^/^Di6, so we must instead produce a Reed-Muller codeword. The code 
TIA4{1,4) is five-dimensional, so we can find a nonzero codeword satisfying any four 
linear conditions. Specifically, we require it to vanish in coordinates 1, 2, ^1 and ^2- 
The resulting codeword c has weight 8 , as do all nonconstant codewords. 

The complementary codeword (1, . . . , 1)— c is also in TlA4(l,4). Its support {^1, . . .,£s} 
contains \,2, k\ and k2 , and it must intersect {^i , . . . ,k%] in exactly four elements. 

Suppose £1 = 1 , ^2 = 2, and £3 and £4 are the other two elements of {£1 , . . . , ^g} \ 
{^1, . . . , ^8}- Now let z = 2~^/^ Y^i=i ^i^ii "^ith ti = t2 = —1 and f,- = —Sj whenever 
li = kj. (The signs and ?4 can be chosen arbitrarily, subject to z having a even 
number of minus signs.) Finally, let w = —2^/^{ti,ei^ + t^e^). As in the previous 
cases, the vectors w,x,y,z span a copy of 2^/^D4 , which completes the proof. □ 
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7 Kissing numbers in dimensions 25 through 3 1 

Before this paper, the best kissing numbers known in dimensions 25 through 31 were 
those shown in Table 3 (achieved by laminated lattices, see Conway and Sloane [15], 
with a correction from Muses [33]). These numbers increase surprisingly slowly from 
each dimension to the next, which made it difficult to believe they could be optimal, but 
they had not been improved since 1982, despite improvements starting in 32 dimensions 
from Edel, Rains and Sloane [22]. 

Dimension Previous record 



25 196656 

26 196848 

27 197142 

28 197736 

29 198506 

30 200046 

31 202692 



Table 3: Previous record kissing numbers 

All of these codes were based on the Leech lattice A24 via a recursive construction in 
which each contains the previous one as a cross section. Recall that A24 is the unique 
densest 24 -dimensional lattice (see Cohn and Kumar [12]) and its kissing configuration 
of 196560 minimal vectors of norm 4 is the unique optimal 24 -dimensional kissing 
arrangement (see Odlyzko and Sloane [36], Levenstein [29] and Bannai and Sloane [2]). 

The old 25 -dimensional kissing number is just 96 more than the value achieved in 
'S?^ , and it is obtained by putting cross polytopes above and below the Leech lattice 
kissing configuration. Specifically, in the coordinates used by Conway and Sloane [19, 
page 133], the Leech lattice contains a copy of 2^/^D24, and to extend its kissing 
configuration to 25 dimensions one simply includes the remaining minimal vectors in 
2^/^D25. Much like the A9 kissing configuration in R*^ (which is constructed in a 
precisely analogous way from ), this code is not even locally jammed: all the points 
not in 2^/^D25 can move. However, the analogy with 9 dimensions is not perfect, and 
the 25 -dimensional version of Proposition 5.1 is not even true. We begin this section 
by showing how to improve the 25 -dimensional kissing number while retaining the 
Leech lattice kissing configuration as a cross section. Once we have shown that, we 
will abandon the cross section constraint and go on to improve all the kissing numbers 
from 25 through 31 dimensions. 

Given a vector v G A24 of norm 6, consider the vectors (2v/3, ±2/\/3) G M^^, 
which lie on the same sphere (of radius 2) as the minimal vectors of the Leech lattice 
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A24 X {0}. If X is any minimal vector in the Leech lattice, then |v — > 4, from 
which it follows that (v,x) < 3. Thus, 

((2v/3, ±2/^/3), (x,0)) <2, 

so all the new vectors stay far enough away from those in the Leech lattice. However, 
they may come too close to each other if we use vectors v that are too close. To avoid 
that problem, we choose a subset R of the vectors of norm 6 in A24 such that for 
v,w £ R with V 7^ w, we have (v, w) < 1 . Then 

((2v/3, ±2/V3), (2w/3, ±2/V3)) < 4/9 + 4/3 < 2, 

as desired. Finally, 

((2v/3, 2/\/3), (2v/3, -2/V3)) =4-6/9- 4/3 < 2. 

Thus, 

{(x, 0) : X e A24, |xp = 4} U {(2v/3, ±2/V3) : v e R} 
achieves a kissing number of 196560 + 2|7?| in M}^ . 

We do not know how to maximize the size of R, but greedy approaches easily yield 
\R\ > 48 via a computer search and thus improve the kissing number. The best we 
have achieved is \R\ = 70, and the vectors are listed in the supplementary infor- 
mation (see Appendix A). If we take into account that all the inner products are in 
{—6, —4, —3, —2, lb 1 , 0} , then linear programming bounds (see Delsarte, Goethals and 
Seidel [20]) prove that \R\ < 280. Thus, our construction of R cannot be improved by 
more than a factor of four, but we do not know how close to optimal it is. 

We do not have a conceptual explanation of the 70-point configuration produced by 
the computer search, but we do have a simple construction that achieves |7?| = 56, 
which is enough to improve the kissing number in R^^ . 

Proposition 7.1 There is an antipodal set R of vectors of norm 6 in A24 such that 
\R\ = 56 and |(x,3')| < 1 for all distinct, non-antipodal x,y ^ R. 

Proof The Leech lattice A24 contains two orthogonal copies of the Coxeter-Todd 
lattice K12 (see Conway and Sloane [19, page 128, equations (129) and (130)]). Thus, 
it suffices to find such a subset of size 28 within Kn- To do so, we will use the 
construction of K12 given in Section 2.5 of Conway and Sloane [16, Section 2.5, 
page 426], which is called A*^^^ in that paper. 

This construction realizes Ku as an Eisenstein lattice, that is, a Z[w]-module with 
u} = e^^'^^ . Let a = 2 + 3uj, so \a\^ = 7. Then the Coxeter-Todd lattice K12 is 
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isometric to L, where 

L = —<.{xi, . . . , xj) G : xi = X2 = ■ ■ ■ = xj (mod a) and ;cy = 1 . 

As before, we use the inner product {x,y) = 2Re(x,3')c, where (z, w)^ = Y^'j=\ ^■i'^j 
is the usual Hermitian inner product on . 

We will make use of only some of the vectors of norm 6 in L, namely those obtained 
by permuting the coordinates of (1, w, w^, 0, 0, 0, 0). Let W be the set of such points. 
Note that \ + uj + up- =0 and a = aoj = aup- = (mod a) , so these vectors are 
indeed in L. 

We begin by finding a six-dimensional regular simplex in Kn (that is, seven vectors 
of norm 6 whose inner products are all — 1 ). To do so, start by choosing seven vectors 
in W such that the supports of any two overlap in exactly one coordinate position. 
Equivalently, label the seven positions by points in the Fano plane ¥^(¥2), and choose 
the seven lines in that plane as the supports of the vectors. 

Suppose two of these vectors overlap in position /, with coordinates oJ and uj'^ . Then 
their inner product in L is 2 Re uJ^'^ , which equals — 1 if and only if j ^ k. Thus, 
we must ensure that each pair of vectors differ in the position where they overlap. 

Equivalently, labeling the incident point-Une pairs in ¥"^(¥2) with 1 , a; or must 
yield a three-coloring of the edges of the bipartite incidence graph between the points 
and lines (the Heawood graph). The existence of such a coloring follows from Konig's 
line coloring theorem (see Lovasz andPlummer [30, Theorem 1.4.18]), which says that 
the chromatic index of a bipartite graph equals its maximal degree. It is not difficult to 
exhibit such a coloring, and hence a six-dimensional simplex in K12, explicitly as the 
rows of the following array: 






1 





OJ 





a;2 





1 








OJ^ 


LO 














OJ 


1 








OJ^ 


OJ 


OJ^ 


1 

















OJ 








OJ^ 





1 


UJ^ 














1 


OJ 








OJ^ 





1 


OJ 






To produce 14 of the desired 28 points in L, we simply take these seven points and their 
antipodes. The remaining 14 can be obtained from them by permuting the coordinates. 
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Specifically, we fix the seventh coordinate and permute the others via the six-cycle 

li-^6i-^4i->3h-^2i-^5i->l 

or its inverse. It is not difficult to check that if we apply either permutation to the 
original simplex in L, then the seven new points all have inner products or ±1 with 
the original points. 

Thus, we have found 28 points in Ku with the desired inner products, and hence 56 
in A24. □ 

In the remainder of this section, we will further improve the records for kissing numbers 
in dimensions 25 through 3 1 . All of our constructions will be quite simple, except 
that they rely on finding a subset S of the minimal vectors in A24 such that no two 
elements of S have inner product greater than 1 . We will give a conceptual argument 
that achieves l^l = 288, which will be enough to improve all the kissing numbers. 
Further optimization using a computer has led to an antipodal subset with l^l = 480, 
which yields the new records shown in Table 1 . The code of size 480 is enumerated in 
the accompanying computer files (see Appendix A). If we take into account that all 
the inner products are in { —4, — 2, it 1 , 0} , then linear programming bounds prove that 
\S\ < [9360/1 IJ = 850. 

Lemma 7.2 There is an antipodal set S of minimal vectors in A24 such that \S\ = 288 
and \ {x,y)\ < 1 for all distinct, non-antipodal x,y £ S. 

The code we construct below is relatively well known, but for completeness we sketch 
the construction here. 

Proof Let MTZ be the Nordstrom-Robinson code, a nonlinear code of size 256 and 
distance 6 in (see Mac Williams and Sloane [31, page 73]). Let W consist of all 
the vectors ( X^,(— l)*'e;)/2 with (si, . . . ,si(i) S MTZ, as well as the vectors ib2e,- 
for 1 < / < 16. The elements of W have norm 4, and is easy to check that distinct 
vectors of W have inner product at most 1 with each other. One can also check that 
the span of W is isomorphic to the Barnes-Wall lattice. Since there is an embedding 
of the Barnes-Wall lattice into the Leech lattice (see Conway and Sloane [15]) we 
obtain the desired subset 5 of A24 as the image of IV. □ 

For the rest of this section, S will denote any set of minimal vectors in A24 such that 
{x,y) < 1 for all distinct x,y £ S. In particular, it may be the set from Lemma 7.2, 
the one of size 480 found by a computer search, or an even larger set. Also, let C 
denote the set of minimal vectors in A24 . 
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Proposition 7.3 The kissing number in dimension 25 is at least 196560 + |5| . 

Proof Fix an angle 9 (to be determined shortly), and let 

T = {(x,0):xeC\S}U {(xcos 9, ±2 sin 9):xeS}. 

In other words, we remove S from the equatorial hyperplane and replace it with two 
copies lying on the parallels at latitude ±9. 

The inner product between a vector (x, 0) and (y cos 9, ±2 sin 9) is {x, y) cos 9 < 2, 
since x / j (and hence {x,y) < 2). 

The inner product between two distinct vectors in the northern hemisphere is 

{(x cos 9, 2 sin 9), (y cos 9, 2 sin 9)) = {x, y) cos^ 9 + 4 sin^ 9 

< cos2 6' + 4sin2 6' 
= 1 + 3 sin^ 9, 

which is at most 2 provided that sin-^ 9 < 1/3 . 

The inner product between two vectors in opposite hemispheres is 



which is at most 2 provided that sin^ 9 > 1/4. Therefore, for any choice of 9 such 
that 1 /4 < sin2 6' < 1/3, the code T of size 196560 + l^l has maximal inner product 



To deal with higher dimensions, we will use disjoint sets 5i, . . . , 551 C C with the 
same property that {x, y) < 1 for all distinct x,y e 5,- . 

Lemma 7.4 Let = S. Then there are disjoint subsets S2,S^, . . . C C such that 
{x,y) < 1 for all distinct x,y G 5,- and 



Furthermore, if S is antipodal, then we can take 52, ^3, . . . to be antipodal as well. 
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{(x cos 9, 2 sin 9), (y cos 9, —2 sin 9)) = {x, y) cos^ 9 — A sir? 9 

< 4 cos^ 9 -A sin^ 9 

< 4 - 8 sin^ 9, 



2 and is thus a kissing configuration. 



□ 
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Of course, Si = if / is sufficiently large. Tiie proof of Lemma 7.4 uses the 
probabilistic method, rather than an explicit construction. This approach is surely not 
optimal, but it is convenient. Furthermore, it naturally leads to a probabilistic algorithm 
that can be used to obtain these sets with the aid of a computer. 

Proof Let G be the automorphism group of the Leech lattice. We will obtain the 
sets S2,St,, . . . by induction as subsets of translates of S by elements of G. First, we 
compute the expected size of the intersection of gS with 5i U • • • U S,- 1 when g is 
chosen uniformly at random from G. Since Si, ... , 5;_i are disjoint, we have 

E 1^5 n (5i u 52 • • • u Si^i)\ = Y.E\8Sn Sj\. 

geG geG 
The terms on the right side are easily computed, via 

E \gSnSj\ = y] Pr lgxC^Sj] 
_ \S\\Sj\ 



196560 

In the last line, we have used the fact that G acts transitively on the 196560 minimal 
vectors; therefore, gx is uniformly distributed among those vectors as g varies, and 

1 

Pr [gx = y] 



geG° 196560 

Now there must be a translate gS whose overlap with U • • • U is at most the 
expectation, namely 

' ' 196560 ■ 

Letting S, = gS \ U'jZi Sj completes the proof. (Note that if S is antipodal, then 
each of these sets will be as well.) □ 

Theorem 7.5 Let Si,S2, . . ■ be as in Lemma 7.4. Suppose we have partitioned a 
kissing configuration in ttie unit sphere S'^~^ C into subsets Ti, . . . ,Tn such that 
for distinct x,y G T,-, we have {x,y) < —1/2. Then the kissing number in dimension 
24 + d is at least 



i=l 
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and it is achieved by the set 



\ {x,0) G X M'' : ;c G C \ |J5,- i U [j {{x^/2j?>,y^/Aj?>) : x G Si,y G T,-}. 

i=\ ^ i=\ 

Note that the minimal angle constraint on implies that |r,| < 3, and |r, | = 3 is 
possible only if T, consists of three points forming an equilateral triangle on a great 
circle. Of course it is pointless to take |r, | = 1 , and to optimize the bound we should 
take iTil > \T2\ > • • • > \Tn\, assuming |5i| > 1521 > • • • > \Sn\- 



Proof First, observe that the points (x-y/2/3,yi/4/3) have norm 

iw' + tM' = i-4 + t-i=4, 



as they should. Within each set { {x^j2fh^y^jAji~) : ;c G 5,, y G T,} , the inner product 
between two distinct points {x^JlJlt^y-yjAf^ and {x! ^Jljl> ^y' ^J~\fi^ is at most 



if J / ' and at most 



2 .4_ 4 . 1 

3 ^ 3 2 



^ • 1 + ^ • 1 

3 ^ ^ 3 ^ 



if ;c / x! . Inner products between two of these sets are also easily dealt with: 
for / 7^ j, the inner product between (A;-y/2/3,yy^4/3) with x ^ Si, y ^ Tj and 
(x' a/273, y a/473) with x' G Sj, y' G Tj is at most 



2.2 + 4.1 

3 ^ ^ 3 2 



2. 



Finally, the inner product between (xy^2/3,3'y^4/3) and (x',0) is at most 2-^/2/3, 
which is strictly less than 2 . □ 



Note that the code from Theorem 7.5 cannot be even locally jammed, provided that 
C / Ur=i '^c because the points in 

|(x,0) G M^^ X : X G C \ U"^') 

i=i ^ 

are not adjacent to any of the points of the code that are not in this set. 

We do not know the best way to apply Theorem 7.5, but we can use it as follows. Given 
any kissing configuration in of size K, we want to partition it into antipodal pairs 
and equilateral triangles (or singletons if necessary). Of course, if it is antipodal we 
can simply partition it into K/2 antipodal pairs, but that will generally not be optimal. 
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Dimension 


Lower bound 




\S\ =288 


\S\ = 480 


25 


196560+ 




196848 


197040 


26 


196560 + 2|5i 1 +2|52| 




197712 


198480 


27 


196560 + 2|5i 1 +2|52| + EL3 




198576 


199912 


28 


196560 + 2 Etil^d 




201156 


204188 


29 


196560 + 2 S^^i|5,-|+Ej=9l 


Si\ 


203430 


207930 


30 


196560 + 2^^^1 




210200 


219008 


31 


196560 + 2 E-=i |5,| + E?i25 


\Si\ 


217588 


230872 



Table 4: New lower bounds on kissing numbers 



If there is an Eisenstein structure, as is the case for A2, D4, E(, and E^, then we 
partition it into regular hexagons using that structure and divide each hexagon into 
two equilateral triangles. For A3 , D5 and £7 , we partition a cross section using its 
Eisenstein structure and then fill the rest with antipodal pairs. This yields the bounds 
shown in Table 4. To compute \Si\, we simply use Lemma 7.4 recursively, taking into 
account that 15,1 must be not just an integer, but also even when Si is antipodal. For 
example, if \S\ = 480 with S antipodal, then \S2\ > 480-(l -480/196560) « 478.83, 
which implies \S2\ =480. 

We do not expect that these bounds are anywhere close to being optimal, and already 
in M^^ they are worse than the previous record: they give 266544, while the record 
is 276032 (l^l = 554 would be required to break the record). However, our bounds 
improve upon all cases from dimension 25 to 3 1 . 

8 Open problems 

We conclude with some open problems. The most basic is whether there are spherical 
codes that span the ambient space and are jammed but not infinitesimally jammed. 
We suspect that there are such codes, but we have not found one. A related question 
is how to test efficiently whether a code is jammed (which would be possible if 
jamming and infinitesimal jamming were equivalent, although one could hope for 
even faster algorithms). Such an algorithm might be based on higher-order variants of 
infinitesimal jamming, but they are far more subtle than one might expect (see Connelly 
and Servatius [14]). 

Our new kissing records in dimensions 25 through 3 1 can presumably be improved, 
and it would be very interesting to know how far these techniques can be pushed, or how 
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to construct much better arrangements. The ratio of the upper bounds from Shtrom [40] 
to the lower bounds proved here grows roughly like 1 .4'^ in 24 + dimensions for 
1 < d < 1 , so there is considerable room for improvement. Note that optimizing the 
sizes of the sets R and S from Section 7 can be viewed as maximizing the size of 
cliques in highly symmetrical graphs. Finding large cliques is NP-hard in general, but 
that does not settle the question of how well one can solve this problem in practice. 

Among kissing problems in low dimensions, dimensions 17 through 23 seem ripe for 
improvement, although it is unlikely that the approach we have used here can be made 
to work. Dimension 16 is particularly interesting, and we would very much like to 
know whether the 4320 minimal vectors of the Barnes-Wall lattice solve the kissing 
problem. They are certainly not the unique solution, because we have found that at 
least one of the 16 -dimensional packings constructed by Conway and Sloane [18] has 
a different kissing configuration of the same size. It seems unlikely that one could 
classify all the possibilities in R'^, but it might be possible to extend the conjecturally 
exhaustive list from Section 4 to higher dimensions, perhaps up through M^^. 

In Section 5 we gave a list of the best kissing configurations known in through 
M^^, but we suspect even more such configurations remain to be found. As mentioned 
in that section, the Ej root system fits into the same framework. In fact its three other 
competitors do too, in the following sense. In each case, one can find a cross polytope 
contained in the configuration (for example, by a randomized, greedy algorithm). Using 
it to define the coordinate system yields simple rational coordinates for each point. In 
the E-j case, these coordinates come from a constant weight code. In the other three 
cases they do not, but the coordinates can be obtained by a systematic mutation of the 
Ej case. Perhaps one could obtain additional kissing configurations in nine through 
twelve dimensions via a similar construction. 

The problem of rigidity naturally generalizes to many other ambient spaces, such as 
projective spaces, but this generalization presents new and interesting features. For 
example, a projective configuration is not determined by its pairwise distance matrix, 
and in complex projective space there are even continuous families of optimal codes 
with exactly the same pairwise distances (see Cohn and Kumar [10, page 129]). 

Even in Euclidean space, which is the most thoroughly studied case so far, there are 
unresolved questions. How can one test whether a periodic packing is jammed? If 
one restricts attention to packings consisting of N translates of a lattice (that is, those 
with particles per unit cell), then one can test for infinitesimal jamming and hence 
jamming (see Donev, Torquato, Stillinger and Connelly [21]). However, the answer 
may depend on A^. For example, the laminated lattice Ag is jammed with = 1 (that 
is, it cannot be deformed as a lattice packing), but not with N = 2 (giving rise to the 
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fluid diamond packings of Conway and Sloane [18]). We know of no bound for how 
large N must be to detect a lack of rigidity. 

Finally, we conjecture that in all sufficiently high dimensions, there exist optimal kissing 
configurations with no contacts whatsoever (that is, no pairs of points with inner product 
1 /2), so they are unjammed in the strongest possible sense. This phenomenon occurs 
in , but we know of no higher-dimensional cases. Part of our motivation for making 
this conjecture is that we know of no large, jammed kissing configurations at all in high 
dimensions. The D„ root system is jammed for « > 3, but it contains only 2n{n — 1) 
points, which is tiny compared with the exponential growth of the kissing number (see, 
for example, Conway and Sloane [19, pages 23-24]). Are there exponentially large 
jammed kissing configurations in high dimensions? Are there even any of greater than 
quadratic size? For example, is the kissing configuration of the Barnes-Wall lattice 
BW2k in M? always jammed? (It has size asymptotic to C • 2*^^*^+^^/^, where C is a 
constant [19, page 24].) 
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A Data files 

As supplementary information for this paper, we have made available sixteen data files 
through the arXiv . org e-print archive, where it is paper number arXiv : 1102 . 5060. 
One can access these files by downloading the source files for the paper. The data files 
can also be downloaded from the web page for this article (doi: 10.2 140/gt.201 1. 15.2235) 

Ten of them describe the kissing configurations enumerated in Section 4. These 
files are each named after the corresponding configuration: 5-40a.txt, 5-40b.txt, 
6-72a.txt through 6-72d.txt, and 7-126a.txt through 7-126d.txt. The first 
line of the file specifies the number N of points, and the second line specifies the 
dimension n of the ambient Euclidean space. The third line consists of n positive 
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integers di, . . . ,dn (separated by spaces), which are the coefficients of the diagonal 
quadratic form used to measure distances. Finally, the remaining A'^ lines each give 
the n coordinates of one of the points in the configuration (again separated by spaces). 
The scaling has been chosen so that all the coordinates will be integers, and the inner 
product between points x and y is defined by 



Equivalently, if we use the standard inner product, then the ith coordinate must be 
scaled by but phrasing it in terms of changing the inner product avoids the need 
to use irrational numbers as coordinates. (There exist different coordinate systems that 
use only rational numbers, even with the standard inner product, but the coordinates 
used here are compatible with the fibering construction from Section 4.) Note that 
within each file, all the vectors have the same norm, but they are not unit vectors. 

Four of the files, namely b9-18.txt, bl0-30.txt, bll-35.txt and bl2-51.txt, 
describe the constant weight binary codes used to build the kissing configurations from 
Table 2 in Section 5. The file bn-N . txt contains all the codes of block length n, size 
N, constant weight 4 and minimal distance 4, up to isomorphism. (We thank Patric 
Ostergard for providing these codes.) Each codeword is given on a line by itself, with 
the binary digits separated by spaces, and there is a blank line between different codes. 

The remaining two files, R . txt and S . txt, describe the 70-point configuration R and 
the 480-point configuration S from Section 7. These files are in a slightly different 
format: they omit the first three lines (N, n and di, . . . ,d„). Instead, each line 
specifies the 24 coordinates of one of the points in the coordinate system used by 
Conway and Sloane [19, page 133, Figure 4.12], but with the irrational factor of 1/ VS 
omitted. 
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